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Abstract 



CJ 

■ Multicast transmission, wherein the same packet must be delivered to multiple receivers, is an important 

'sj" ' aspect of sensor and tactical networks and has several distinctive traits as opposed to more commonly studied 

(N : 

unicast networks. Specially, these include (i) identical packets must be delivered successfully to several nodes. 



(ii) outage at any receiver requires the packet to be retransmitted at least to that receiver, and (iii) the multicast 
rate is dominated by the receiver with the weakest link in order to minimize outage and retransmission. A first 
contribution of this paper is the development of a tractable multicast model and throughput metric that captures 



j. I each of these key traits in a multicast wireless network. We utilize a Poisson cluster process (PCP) consisting 

> 



of a distinct Poisson point process (PPP) for the transmitters and receivers, and then define the multicast 



jn ■ transmission capacity (MTC) as the maximum achievable multicast rate per transmission attempt times the 

' maximum intensity of multicast clusters under decoding delay and multicast outage constraints. A multicast 



cluster is a contiguous area over which a packet is multicasted, and to reduce outage it can be tessellated into 



(N 
O . 

. V smaller regions of multicast. The second contribution of the paper is the analysis of several key aspects of 



this model, for which we develop the following main result. Assuming t/v transmission attempts are allowed 
for each tessellated region in a multicast cluster, we show that the MTC is O {pk^ log(fc)w^) where p, x and 



^1 y are functions of t and v depending on the network size and intensity, and k is the average number of the 

intended receivers in a cluster We derive {p, x, y} for a number of regimes of interest, and also show that an 
appropriate number of retransmissions can significantly enhance the MTC. 
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I. Introduction 

Multicast refers to the scenario whereby a transmitter needs to send a packet to multiple receivers. 
In a wireless network, this creates a two-edged sword. On one hand, the broadcast nature of wireless 
transmission assists multicast; but roughly uncorrected outage probabilities at each receiver (due to 
spatially distinct fading and interference) require retransmissions that cause interference and waste. 
Multicast is an important aspect of sensor and tactical networks, and increasingly in commercial 
networks where streaming is supported. However, the literature on multicast is minuscule compared 
to unicast - whereby nodes are paired into sources and destinations - and even basic modeling issues 
such as outage, capacity/throughput definitions, and retransmissions are not widely agreed upon. In 
this work, we attempt to investigate the fundamental throughput limits of multicast transmission and 
we develop a metric based on spatial outage capacity which we term multicast transmission capacity 
(MTC). 

In order to characterize the MTC in a wireless network we propose a multicast network model 
in which each transmitter has an intended multicast region (called a cluster) where all the intended 
receivers are uniformly and independently scattered, and hence a Poisson cluster process can be 
reasonably used to model the transmit-receiver location statistics. The active transmitters are modeled 
as a stationary Poisson point process (PPP) and their associated receiver nodes in the cluster are also 
a stationary PPP, as shown in Fig. [TJ In other words, each cluster is randomly located in the network 
and comprises a multicast session. This paper will develop interference and outage expressions for 
this model, and analyze some important cases of the network model and design space, including the 
network intensity, size and the effect of retransmissions. 

A. Motivation and Related Work 

The majority of the existing works on network capacity are focused on the unicast scenario 
and built upon the protocol and physical network models proposed in Generally speaking, the 
unicast network capacity is the maximum number of point-to-point communication links that can be 
simultaneously supported in the network under some transmission constraints. For example, transport 
capacity built upon the protocol model has a geometric constraint on transmitter-receiver pairs. Due 
to the difficulty in coordinating the transmission constraints between multiple receivers, the unicast 
capacity, in general, is not readily extended to the multicast capacity. Even the definition of multicast 
capacity is not widely agreed upon. 
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Some previous works, such as [|2ll- [fTT]| . have made significant progress in studying the multicast 
or broadcast capacitj^. For example, in [|2l the protocol model is used where source nodes and their 
multicast destinations are randomly chosen. The multicast capacity is defined as the sum rate of 
all multicast flows and it is obtained as a function of the number of multicast sources. In [4], the 
multicast capacity under the protocol model is defined as the transmission rate summed over all 
of the multicast traffic flows in the network. Its scaling characterization is obtained by the number 
of receivers in each multicast session. Reference [|9]| showed that the broadcast capacity under the 
protocol model does not change by more than a constant factor when the number of nodes, the radio 
range or the area of the network is changed. In fTOl, the physical model and a stationary PPP of the 
nodes in the network are considered. It showed that the broadcast capacity is a constant factor of the 
computed upper bound when the number of nodes goes to infinity under a constant node intensity. 
The multicast and broadcast capacities in the above works are defined behind the main concept that 
the transmitted information should be received by all of its intended receivers. However, they are not 
investigated from the multi-receiver outage point of view and thus their scaling results cannot provide 
us retransmission guidelines for capacity enhancement. 

The multicast capacity problem in this paper is studied from an outage perspective, which is a 
departure from previous work. For example, the prior work built on the protocol model in [1] does not 
consider channel impairments such as fading and path loss, and thus, all transmissions are successful 
once certain geometric constraints are satisfied. Also, outage and the resulting retransmissions are 
not considered. This is somewhat unrealistic. The main issues we would like to clarify are (i) how 
many simultaneous multicast sessions can (and should) coexist when all receivers in a multicast 
session need to receive the packet from their transmitter, and (ii) when are retransmissions beneficial 
or detrimental to the multicast capacity? Hence, we introduce the MTC in a planar network, which 
is defined as the maximum achievable multicast rate per transmission attempt times the maximum 
number of the coexisting clusters in the network per unit area, subject to decoding delay and multicast 
outage constraints. This is a logical evolution of the transmission capacity framework originated in 
[fT2l to the topic of multicast. The decoding delay constraint here means a transmitter can multicast 
a packet to all of its intended receivers up to r G N+ transmission attempts and multicast outage 



'Broadcast means only one transmitter would like to transmit to all of other nodes in the network, whereas multicast means transmitters 
transmit to a certain number of nodes in the network. Hence, broadcast capacity can be viewed as a special case of multicast capacity. 
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happens when any of the intended receivers in a cluster does not receive the information multicasted 
by their transmitter during those r attempts. 

B. Main Contributions 

Since unicast outage cannot be directly applied to the outage scenario of multiple receivers, our 
first contribution is introducing MTC with multicast outage. To the best of our knowledge, this is the 
first work to study multicast capacity under a multicast outage constraint. According to the proposed 
Poisson cluster model for multicast, we also propose some new cluster-based definitions for the 
largeness and denseness of a network in order to characterize the scaling behaviors of the MTC under 
different conditions. Referring to Table IJ for the notation of main variables, we have shown that the 
scaling of both single-hop and multihop MTC can be expressed in a general form of Q(pk^ log(A;)t>^) 
where p, x and y are given in Table |ll] for different network conditions, k is the average number of 
the intended receivers in a cluster, and v is the number of the tessellated regions of equal area in a 
clustej^. From Table HIl we know retransmissions have a significant effect on the MTC and certain 
number of retransmissions could enhance it. We found that the decoding delay r can be viewed as 
a resource which should be allocated properly in every tesselated region to maximize the MTC. In 
addition, we also show that the MTC scaling holds for the various parameters of Nakagami fading. 

We characterize three approaches to improving the MTC - interference-suppression, interference- 
avoidance and area- shrinking (i.e. reducing the cluster size). The area- shrinking method is able to 
provide the best capacity gain among the three since the multicast outage probability is reduced due 
to fewer receivers in a cluster. However, shrinking the area of a cluster means some of the intended 
receivers may have to be excluded, which is not welcome. Thus the multihop multicast method is 
proposed to improve the MTC without shrinking the cluster. Our main result shows that if the clusters 
are appropriately tessellated multihop multicast can significantly improve the MTC compared to its 
single-hop counterpart. This is because an appropriate number of retransmissions largely reduces the 
multicast outage probability and thus results in the increase of the maximum contention intensity that 
compensates the loss of spectral efficiency due to retransmissions. 

^Throughout this paper, standard asymptotic notation will be used: O(-), and 0(-) correspond to asymptotic upper, lower, and 
tight bounds, respectively. 
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C. Paper Organization 

In Section [III the network model and its assumptions are described and some preliminaries for 
the following analysis are provided here. The main results of multicast transmission capacity for the 
case of single-hop multicast are presented in Section |llll whereas Section |IV] has the main results of 
multicast transmission capacity with multihop multicast. Finally, we conclude our work in Section |Vl 

TABLE I 

Notation of Main Variables, Processes and Functions 



Symbol 


Definition 


$' ($■", $'^) 


PPP of transmitters (receivers, connected receivers) 


e 


upper bound of multicast outage probability (<C 1) 


At (Ar, Ac) 


intensity of <I>' , ^'^) 


At 


maximum contention intensity of At 


A. 


first-order Taylor expansion of At (Aj ~ At for small e) 




multicast outage event depending on At 


T 


decoding delay constraint, positive integer 


V 


number of the tessellated regions in a cluster (v < r) 


s 


radius of a multicast cluster (s > 1) 




Lebesgue measure of a bounded set A 




SIR threshold 


a 


path loss exponent (a > 2) 


b 


multicast transmission rate (bps/Hz) 


a 


multicast transmission capacity (bps/Hz/m'^) 


k 


average number of the intended receivers in a cluster, vrs^Ar 


/„(.), Fh(-),F«(-) 


PDF, CDF, CCDF of channel fading gain H 



TABLE II 

Summary of Main Results on Multicast Transmission Capacity 



Multicast Transmission Capacity d 




Maximum Contention Intensity A^ 




Network Condition 


Dense 


Large 


Large Dense 


X 


V 
T 


-(1 + 7) 


-(^) 


P 




y 


7 + 1 
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II. Network Model and Preliminaries 

In most of the existing literature, multicast capacity is defined based on the sum of the supportable 
rates of all multicast sessions. Similarly, the MTC here is characterized by the sum of the multicast 
rates of the multicast sessions not in (multicast) outage. The multicast rate is affected by the number 
of the receivers as well as the interference from other multicast sessions. The multicast network model 
is developed with the principle that it should capture the key traits of a multicast network (for example 
spatial reuse, retransmissions, broadcast packets, etc.), while being as tractable as possible. 

A. Clustered Network Model for Multicast Transmission 

In the network, each transmitter has a multicast cluster of equal area and its receive nodes in the 
cluster suffer aggregate interference from a Poisson field of transmitters. Specifically, we assume that 
the network is operating a slotted ALOHA protocol and the distribution of the transmitting nodes in 
the network is a stationary Poisson point process (PPP) of intensity At. As shown in Fig. [H any 
transmitter Xj G has its own intended multicast cluster TZj where all of its intended receivers are 
uniformly and independently distributed and they also form a stationary PPP $• of intensity Ar. Note 
that each cluster could contain other transmitters and unintended receivers in addition to its own 
transmitter and intended receivers. 

Accordingly, the multicast transmission sessions in the network follow a the Poisson cluster process 
(PCP) Zi = $^UXj, i.e., each transmitter is a parent node associated with a cluster of receive daughter 
nodes. The cluster processes [Zi] corresponding to different transmitters {Xi} are assumed to be 
independent so that the superposition of all clusters yields the resulting cluster process $ = IJx e*' ^« 
of intensity \ = k\\ where k = Trs^Ar is the average number of the intended receivers in each cluster 
assuming all {T^j, Vz G N} have the same radius s > 1. The distribution of the intended receiver nodes 
in each cluster is modeled as a marked PPP denoted by $^ = {{Yij,Hij) : Yij G TZi,j G N}, where 
Hij is the fading channel gain between transmitter Xi and its intended receiver Yij. Similarly, the 
distribution of transmitters in the network is also a marked PPP, i.e., = {{Xi, {Hij}),i,j G N+} 
where Hij denotes the fading channel gain between transmitter Xi and the receiver Yqj located in 
cluster T^-ol^ All fading channel gains are i.i.d. with probability density function (PDF) 

^ Since all of the following analysis is based on the nodes in the reference cluster TZo, the subscript of some variables will not be 
explicitly indicated if there is no ambiguity. So Y, and Hj in TZo actually stand for Yoj and Hoj, respectively. 
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Without loss of generality, the MTC can be evaluated in the reference cluster TZq whose transmitter 
Xq is located at the origin. We condition on this typical transmitter Xq resulting in what is known the 
Palm distribution for transmitting nodes in the two-dimensional Euclidean space [fT3l . It follows by 
Slivnyak's theorem [fT3l that this conditional distribution also corresponds to a homogenous PPP with 
the same intensity and an additional point at the origin. The signal propagation in space is assumed 
to undergo path loss and fading. The path loss model between two nodes X and Y used in this paper 

x-r if x-r > 1 
ii\x-Y\)^r (1) 

I 0, else, 

where |X — F| denotes the Euclidean distance between nodes X and Y, and a > 2 is the path 
loss exponenj^. The reason of using the model in ([T]) is because the model | ■ |~" does not behave 
well in the near field of each transmitter and it thus leads to an unbounded mean of the shot noise 
process. This model is similar to the idea of the bounded propagation model proposed in [fT6l [fTTll . 
The Nakagami-m fading model is adopted in this paper because it covers several different fading 
models, such as Rayleigh (for m = 1), Rician fading with parameter K (for m = (K + 1)^/ {2K + 1)) 
and no fading (for m — t- oo), as well as intermediate fading distributions. It is of unit mean and 
variance and given by 

fHih) = :^h"^-'expi-mh), (2) 
1 [m) 

where m is a positive integer and T{m) = z'^~^e~^ dz is the Gamma function. 

Each receiver is able to successfully receive its desired information if its SIR is greater or equal 
to the target threshold (3. That is, receiver node Yj is "connected" to the typical transmitter if 

SIR ^ J^l^iM > /3. (3) 

Note that SIR depends on At, i.e. SIR = SIR(At). All the transmitters are assumed to use the same 
transmit power, the network is interference-limited, and Ij is the aggregate interference at receive 
node Yj and a sum over the marked point processes. Namely, 

/,= Yl H,,ii\X,-Y,\), (4) 
x,e'S>\{Xo} 



''in a planar network, a is greater than 2 in order to have bounded interference, i.e. Ij < oo almost surely if a > 2 1141 . 1151 
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which is a Poisson shot noise process, and Hij is the fading channel gain from transmitter Xj to 
receiver Yj in TZq. Since is stationary, according to Slivnyak's theorem the statistics of signal 
reception seen by receiver Yj is the same as that seen by any other receivers in the same cluster. Thus 
Ij can be evaluated at the origin, i.e., dH) can be rewritten as Jq = J2xie'S>t\{Xo} Hii{\Xi\), where Hi 
is the fading channel gain between transmitter Xj and the origin. 

Suppose the decoding delay is up to the lapse of r transmission attempts for a transmitter. The 
connected receiver process for the t-th transmission is denoted by 

^it) = {(Y^H^it)) e $[, : H^miYjl) > m, (5) 

where {Hj(t)} are i.i.d. for all t E [1, ■ ■ ■ ,t]. Also, let $^(r) be the connected receiver process at 
the rth attempt, i.e., it is the set of all intended receivers in a cluster connected by their transmitter 
during the decoding delay, and thus it can be written as $*^(r) = IJt=i other words, the 
connected receiver process can be described by a. filtration procesjfi 

B. Multicast Transmission Outage 

The transmission capacity of an ad hoc network introduced in [ [T2l is defined based on point-to-point 
transmission with an outage probability constraint e G (0, 1), and is given by 

c, = 6At(l-e), (6) 

where b is the constant transmission rate a communication link can support (for example, about 
log2(l + (3)), and At is the maximum contention intensity subject to an outage probability target e. 
However, Q cannot be directly applied to multicast because the multicast rate would be affected by 
Ar and At, and the outage of a multicast transmission is not point-to-point but point-to-multipoint. 
How to declare an outage event for a transmitter multicasting information in the previous multicast 
transmission model is a key issue. 

Since no desired receiver can be assumed to be dispensable, a reasonable way to define multicast 
outage is when any of the intended receivers of a transmitter does not receive a multicasted packet 
during a period of time up to the decoding delay. That is, after all the allowed retransmissions have 
been used, if one of the desired receivers in the cluster has not decoded the packet, we declare an 

'a filtration process means $'^(1) C $^^(2) ■ ■ ■ C <l?'^(r), and for any set A C TZo, A{^'^{t)) .4("l>o) almost surely as r — s- cx) 
where A{^) denotes the random number of point process <1> enclosed in set A. 
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outage for this cluster. Thus, a mukicast outage event of each muMcast cluster can be described as 
£{\i) = {$[) \ $*^(r) 7^ 0} because £ depends on At. The probability of £{Xi) can be characterized 
by the intensity of the connected receivers during the lapse of r attempts as follows: 

P[£(A,)] ^ 1-P[{<l>^\$^(r)} = 0] 

= l-exp|-^ (Ar-Ac(r,r))Mdr)|, (7) 

where //(■) denotes a Lebesgue measure and Xc{Y, r) is the intensity of $'^(r) at node Y. Using ^ to 
find multicast outage probability can be interpreted as finding the void probability of a "disconnected" 
PPP in a cluster. Since all of the intended receivers are uniformly distributed in T^o, © can be rewritten 
as 

P[£(At)] = l-exp{-7rs2(Ar-EK[Ac(i?,r)])} <e, (8) 

where Re [0, s] is a random variable whose PDF is //j(r) = ^. The outage probability in ([8]) cannot 
exceed its designated upper bound e which is assumed to be a small value throughout this paper. 

Remark 1. In Section I7//-AI we will show that $'^(r) is a nonhomogeneous PPP and the average of its 
intensity E/j[Ac(-R, t)] can be found. Thus, the multicast outage probability in ([8]) can be calculated. 
Also, E/j[Ac(-R, t)] is a monotonically decreasing function of Xt so that the maximum Af is attained 
when Kr[Xc{R, t)] reduces to its lower bound Xr + ^"^^^^2'^^ th<^t is obtained by solving ([8]). 

Remark 2. The multicast outage probability in ([8]) will accurately approximate the unicast outage 
probability for a small e if there are no retransmissions (r = 1) and only one receiver is in a cluster. 
For a unicast scenario in each cluster, we have ns'^Xr = 1 and E^[Ac(-R, 1)] = A,-P[SIR(Af) > (3]. 
Thus, dH]) becomes 

F[S{Xt)] = 1 - exp(-P[SIR(Af) < /?]) = P[SIR(Af) <f3] + 0{e^) 

since P[SIR(Af) < f3] < e and e"^ = 1 — e + 0(e^). This shows that the point-to-point outage scenario 
is covered by our model, i.e. multicast outage is a generalization of point-to-point outage. 

C. Definitions of Multicast Transmission Capacity, Largeness and Denseness of Networks 

Since the multicast outage probability is upper bounded by a small e, the maximum contention 
intensity At is a function of e and its Taylor expansion for e gives 

A,(e) ^ sup{At > : P[£:(At)] < e} = A, + 0{e^), (9) 
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where is the Taylor expansion of At(e) without the second and higher order terms of e. Since e is 
small, Ae ~ At(e) and thus, for simplicity, we will focus the analysis on A^ in the following. 

Definition 1 (Multicast Transmission Capacity). The multicast transmission capacity with the 
multicast outage probability defined in ^ for small e is defined as 

a = -&A,(l-e), (10) 
r 

where \ is the first order approximation of Af(e) as indicated in h is the maximum achievable 
multicast rate on average for every cluster and it is not a constant in general. 

Multicast transmission capacity gives the number of successful multicast clusters with the 
maximum achievable multicast rate, that can coexist per unit area subject to decoding delay and 
multicast outage constraints. In other words, it is the area spectral efficiency of cluster-based multicast 
transmission. The following definitions of largeness and denseness of a network will be needed to 
acquire the scaling characterizations of the MTCs in the subsequent analysis. They are defined based 
on the circumstance that the average number of the intended receivers in a cluster is sufficiently large, 
i.e., k = TTs^Ar > 1. 

Definition 2 (Denseness and Largeness of a network with a PCP). (a) We say a network is "large" 
if the area tcs"^ of a cluster in the network is sufficiently large such that for a fixed intended receiver 
intensity A^ we have k ^ 1. (b) If the intended receiver intensity is sufficiently large such that for 
fixed area vrs^ we have A; ^ 1, then such a network is called "dense", (c) A "large dense" network, 
it means that clusters in a network have a sufficiently large size as well as intended receiver intensity; 
namely, \r oc tts^ and thus k ^ 1. 

Here we should point out that Definition [2] may not be consistent with some popular node- 
based unicast definitions in prior literature. For example, a dense network usually means it is dense 
everywhere (i.e., uniformly dense); however, our denseness definition could involve the case of local 
denseness if the receiver intensity in a single cluster is sufficiently large whereas the cluster intensity 
is small. 

D. Multicast Transmission Methods - Single-hop and Multihop 

The network model with a PCP for multicast introduced in the previous subsection implicitly 
assumes that the parent node of each cluster (see Fig. d]) is the sole transmitter. This is the case of 
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single-hop multicast. Hence, for single-hop multicast, multicast outage probability only depends on 
the channel conditions between the parent node and its intended receiver nodes. When the size of 
clusters is large, the path loss of transmitted signals and the average number of the intended receivers 
for each parent node are both increased significantly. The MTC in this case will correspondingly 
decrease, so single-hop multicast is not an efficient means of disseminating information for a large 
and/or dense network. 

To alleviate this drawback, we propose a multihop multicast approach. The idea is to allow 
retransmissions in a cluster by randomly selected receivers that have successfully received the packet 
already. Each of the selected receivers has its own small local multicast region, and the whole cluster 
is covered by the combination of small multicast regions. Note that only one selected receiver is 
allowed to transmit for each time slot in order to make all transmitters in each time slot still form a 
PPP. The detailed algorithm and modeling assumptions will be presented in Section IIV-AI 

III. Multicast Transmission Capacity with Single-Hop Multicast 

In this section, we study the MTC when transmitters are multicasting to all of their intended 
receivers in a single-hop fashion. First we have to find the multicast outage probability defined in 
([8]) and thus we need to study the intensity of the receiver-connected process during the lapse of r 
attempts. Then the maximum contention intensity which characterizes the single-hop MTC can be 
found based on the intensity of the connected receivers in a cluster. 

A. The Receiver-Connected Process 

During the allowed transmission r attempts, the intended receivers in IZq connected by the transmit- 
ter Xq form a receiver-connected process whose intensity is the necessary information to estimate the 
multicast outage probability. Since $° is a filtration process and upper bounded by (as explained 
in Section III-AI) . the connected receiver intensity in IZq is an increasing function of r as shown in 
the following lemma. 

Lemma 1. Consider the stationary PPP $q of intensity in the reference cluster TZq. If a transmitter 
is allowed to transmit a packet up to r times, then $'^(r) is a nonhomogeneous thinning PPP and the 
intensities of ^'-^{r) for different fading models are shown as follows. For Rayleigh fading, we have 

Ac(r, r) = A,(l - {1 - exp [-7rAi(/3r^ oo)Xt]V) , (11) 
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and Ai(-, •) is defined in Proposition [7] in Appendix \I-A\ For Nakagami-m fading with m > 1, we 
have 

Ac(r,r) = A,|l - [1 - ^("-^^(m/? r")] , r G [l,s] (12) 

where 

vi/M(0) A y__lJP exp|_^AfAi(0,oo)-log0}. (13) 

m! a0™ 

Proo/.- See Appendix III-A[ □ 
The connection intensity Ac for Rayleigh fading in (fTTI) reveals an interesting implication. The term 
1 — exp[— 7rAi(/3r", oo)At] can be interpreted as the probability that there is at least one interferer in 



the circular area of radius A/Ai(/3r", oo). This circular area can be called the dominating interferer 
area centered ar r because any single interferer within this area can cause an outage at the receiver 
located at r.The effect of retransmission can be said either to make this probability reduce by r- 
fold or to enlarge the circular area. In a dense network with At > Ar, for example, the term (1 — 

exp[— 7rAi(/3r", oo)At])'^ ^ 1 — r exp[— 7rAi(/3r", oo)At] since 7rAi(/3r", oo)At ^ 1. So the radius of 



the dominating interferer coverage is approximately increased by A/ln(r)-fold if r > 2 in this case. 

In addition, as shown in (nT)) -(fT3l). we know that a closed-form expression of the average connection 
intensity E[Ac] in terms of At in a general network is difficult to find. However, an upper bound on 
E[Ac] can be found in the following lemma, and thus in the special case of clusters with many intended 
receivers, ([TT]) and (fT2]) can be simplified to allow a nearly closed-form solution of At (see Section 
HlFB]) . 

Lemma 2. The connection intensity Ac of a non-homogeneous PPP $'^(r) in a cluster for Rayleigh 
fading is shown in (fTT)) and its expression for Nakagami-m fading with m > 1 is given by (fT2)) . The 
upper bounds on Xcfor different fading models are given in the following: 

Er[Xc{R,t)] < Xr[l- il-ER[exp{-nXtAi{l3R'',oo))]y], for Rayleigh fading (14) 
ER[XciR, ^)] < A, [1 - (1 - Ejj[^('"-^)(m/3/2")])^] , for Nakagami-m fading, m > 1. (15) 

Proof: According to the Holder inequality, for two real-valued random variables Di and D2, 
(E[DiD2]y < E[D^^]E[Dl]P/i, where p, g > and 1 /p+l/q = 1. Consider the number of transmission 
attempts r > 1 and = 1, and let p = r and q = Then it follows that {E[Di]y < E[Dl]. For 
Rayleigh fading, taking average on the both sides of (fTTI) and using this property (E[Di])^ < E[£)[] 
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by letting Di = 1 — exp(— 7rAtAi(/3i?", oo)), (fT4l) follows. Similarly, (fT5l) can be shown in the same 
way. □ 

B. Single-hop Multicast Transmission Capacity 

Now we characterize the MTC in a network when a transmitter directly multicasts its intended 
receivers in a cluster. 

Theorem 1. Suppose the multicast outage probability given in ([8]) is upper bounded by small e and the 
maximum decoding delay is r transmission attempts. If the average number of the intended receivers 
in a cluster with radius s is k and k > -f^, then the maximum contention intensity is 

A. = ^^ = ef^^VV (16) 

where p = ^ \J e(r + 1), /3 is the SIR threshold for successfully decoding and rj is a constant 
(depending on m, (5 and a). 

Proof: See Appendix III-Bi □ 

Remark 3. The scaling function ©(•) ofX^ in (fT6l) only contains the "controllable" network param- 
eters such as s, e, r and (3, which means their values are adjustable if needed. Constant rj contains 
the parameter m of Nakagami fading, which is a channel characteristic and usually uncontrollable 
and thus rj is not left in ©(■). 



If a unicast planar network without retransmission is considered (i.e. k = t = 1), in (1161) will 
reduce to the previous results discovered, i.e., = 6 ^^r^y^j. In [fT2|. for example, the maximum 
contention intensities of FH-CDMA and DS-CDMA are 6 (^-JpjT^^ and 6 (^f^^) respectively, 
where M is the channel number of FH-CDMA and the spreading factor of DS-CDMA. It is easy 
to check that these two results coincide with ours here by considering |y for FH-CDMA and for 
DS-CDMA. In addition, the longest transmission distance in a cluster is s and we know = 6(s^^) 
and so is the network capacity, which also coincides with the results in [f T2| [fT4ll . 

The result in (fT6l) only indicates how much the maximum intensity of transmitters can be supported 
in a network under the decoding delay and multicast outage constraints. Having the maximum 
contention intensity only is unable to tell us how much its corresponding MTC should be since 
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the mukicast rate b is also affected by the maximum contention intensity. Considering a capacity- 
approaching code is used, the maximum achievable multicast rate b that is acceptable for all intended 
receivers is the following ergodic channel capacity evaluated at the boundary of a cluster: 



E 



loe 



(17) 



where -f^max = ^^^te[i,- ,t] H{t) and H{t) is the fading channel gain for the t-th transmission between 
typical transmitter Xq and a receiver located on the cluster boundary. Although there may be no 
receivers on the boundary of a cluster, multicast rate b should be considered from a worse case point 
of view because it needs to be acceptable for all intended receivers in any locations within a cluster. 
The bounds on b are given in the following lemma. 

Lemma 3. There exists a 5 G (0, 1) such that the bounds on the multicast rate b in (fTTl) can be given 
by 



log ( 1 + ) < & < log ( 1 + ] + 0(1). (18) 



Proof: See Appendix Ill-Ci □ 



Remark 4. The bounds on the multicast rate in (fTSl) are not affected by channel fading because the 
fading effect has been averaged out. Lemma\3\suggests that b is significantly reduced by the aggregate 
interference from the transmitters in the cluster and the transmitters out of the cluster only reduce it 
by at most a constant. 

Scaling Law of Single-hop MTC. According to Theorem \T\ and Lemma [3l we found that the 
multicast rate b is 0(log(A;)/r) for any network conditions if is achieved and k is sufficiently 
large. The MTCs in a network without receiver cooperation can be concluded as follows, (i) For 
a dense network, Ae = O (^^j since tts^ is fixed and A; ^ 1. By the MTC definition, we know 
C, = Q (^^^). (ii) If the network is large, then A, = 6 (^^) and thus C, = Q (^^^^ 
(iii) For a large dense network, Ar = 6(s^) and A; ^ 1. So A^ is 6 ^pr^ A;^(^) j and thus 
9 ^pA;^(^) log(A;)j. In summary, the MTC here can be expressed in a general form as follows: 

c, = eipk^\ogik)), (19) 

where x has been given in Table Ull Note again that the scaling function 6(-) is applied only to the 
main controllable variables for a specific network. For example, the main controllable parameters for 
a dense network are Ar (or k) and r whereas vrs^ is a constant. 
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A simulation example of the MTCs for a large network with Rayleigh fading is presented in Fig. 
[2l One can see that the MTCs decrease when s increases (i.e. k increases), and slightly increasing 
the radius of a cluster can significantly reduces the MTC. The decoding delay constraint r has a 
significant effect on the MTCs as well. Since retransmissions decrease outage probability as well as 
increase interference, it can be observed from (flOl) that there exists an optimal tradeoff between r 
and Ae- This can be observed in Fig. [21 for example, the MTC of r = 3 is larger than that of r = 10, 
which indicates a few retransmissions indeed increase the MTC and too many retransmissions are, 
on the contrary, detrimental to it. Fig. [3] presents the MTC results for a large dense network with 
Rayleigh fading. If we compare Fig. [3] with Fig. [21 we can see that denseness does not have a serious 
impact on MTC as largeness. This observation coincides with the scaling law stated in above. Thus, 
path loss is the main key issue of limiting the MTC, which enlightens us the idea of using multihop 
multicast to deliver packets with a less path loss (see Section [TVl). 

IV. Multicast Transmission Capacity with Multihop Multicast 

The MTCs for single-hop multicast investigated in Section [Hlj scale like 0(s^^) if other network 
parameters are fixed. Therefore, the single-hop MTC increases when its multicast cluster is shrunk. 
However, shrinking the cluster is not welcome if the packets must be transmitted over the same 
coverage. So here we would like to know if there is another method to increase the MTC without 
shrinking the cluster. From previous results, for sufficiently large k the scaling of the single-hop MTC 
can be written as 

c.^e(^^^), (20) 

where p = ^{/eij + 1) and x = —K So ([201 ) suggests three approaches to increasing the MTC: 
interference-avoidance, interference-suppression and area- shrinking methods. The capacity gain due 
to interference avoidance can be acquired by removing co-channel interferers such that SIR(At) is 
improved by reducing At. This is the context when each transmitter independently selects its own 
transmission channel from several available channels (see the case of FH-CDMA in |[T2l '). Interference 
avoidance does not affect any network parameters except the multicast rate b. The interference- 
suppression capacity gain can be obtained by signal processing techniques to increase the SIR so that 
multicast rate b is increased (see the case of DS-CDMA in |[T2l '). In addition, suppressing interference 
is equivalent to relaxing the SIR threshold /3 and thus pre-constant r] depending on f3 is increased. 
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Of course, the area- shrinking capacity gain is attained by multicasting in a smaller region instead of 
the whole cluster. Shrinking a cluster only leads to a decrease in the cluster radius s and thus the 
average number of the intended receivers k becomes smaller. 

Let {ga,gs,gv} > l be the interference-avoidance, interference-suppression and area- shrinking gain 
parameters, respectively. Then respectively replacing s^, (3 and by s'^/g^, (3/ Qs and \e/ ga in (l20l) . 
we obtain 

The three gain parameters in (|2T]) indicate which method is able to contribute more to the MTC in 
each context. The interference- avoidance method is superior to interference-suppression when Qa and 
Qs are equal because ga/gV"' > 1 due to a > 2. This point has been shown for FH-CDMA and 
DS-CDMA in [fT2|. If the cluster is shrunk to a smaller region of area ns'^/g^ and all the three gain 
parameters are equal, then shrinking is the best way to improve the MTC. According to (|2T1) . we can 
conjecture that the MTC could be increased if the cluster is tessellated into several smaller regions 
and a packet is allowed to be multicasted in each of them up to some times under the condition that 
the decoding delay and multicast outage constraints both have to be satisfied. This conjecture will be 
verified later in the following subsection. 

A. Multicast over Multihop 

Suppose the cluster IZq is tessellated into a certain number of smaller multicast regions of equal 
area. A packet is multicasted the same number of times (called a multicast time slot) in each tessellated 
region. So the packet is delivered slot by slot from the central typical transmitter to those regions 
in a certain order. Note that the number of the tesselated regions cannot exceed the decoding delay 
constraint r, i.e., the packet is delivered to all its intended receivers at most r — 1 hops. As shown 
in Fig. m for example, each cluster consists of 6 smaller tessellated regions and the decoding delay 
constraint for each region is 2 if r = 12. 

Delivering a packet by the multihop multicast method proceeds as follows. In the first time slot, the 
typical transmitter aims at multicasting its own region. In the second time slot, the typical transmitter 
randomly selects a receiver in the neighboring regions that successfully received the packet in the 
previous time slot, then the selected receiver becomes the transmitter for the next time slot. The 
multihop multicast method proceeds slot by slot in this way until all the tessellated regions are 
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visited by the packet. Note that in each time slot the transmitter is asked to ensure all the receivers in 
its region and at least one receiver in the neighboring regions should receive the packet; otherwise, 
there is an outage. In addition, although there might be more than one receivers in the neighboring 
regions which successfully receive the packet; however, only one of them is selected to multicast in 
the next time slot in order to satisfy the assumption of the PPP of the transmitters (see Proposition [2] in 
Appendix II-BI for the dual property between PCP and PPP). The scaling realizations of the maximum 
contention intensity for the above multihop multicast method are shown in the following theorem. 

Theorem 2. Suppose cluster TZq is tessellated into v smaller multicast regions of equal area and 
each packet is allowed to be multicasted in each region at most r/v times. If the average number 
of the intended receivers in a cluster k > py^, then the following scaling characterization of the 
maximum contention intensity \ is achieved with high probability: 



where x = y = ^ + 1, p = -^{e{r/v + l))r and rj is a constant depending on m, (5 and a. 



By comparing (fT6l) and (l22l) . it is hard to see if the multihop multicast method achieves a larger 
Ae- Nevertheless, in the following subsection we will show that and its corresponding MTC indeed 
achieve a larger value by multihop multicast. The simulation result of the MTC with multihop multicast 
in a large network with Rayleigh fading and Ar = 0.2 is presented in Fig. |5l and it is easy to observe 
that tessellating the cluster into a certain number of regions improves the MTC whereas too many 
tessellations degrades it. Fig. [6] shows the simulation of the MTC achieved by multihop multicast 
in a large dense network with Rayleigh fading. If we compare Fig. |6] with Fig. |5l we can find that 
the curves in these two figures are not very much different. This point is quite different from the 
case of single-hop multicast (see Figs. [2] and [3]), and therefore, it reveals that multihop multicast is 
able to efficiently alleviate the impact on MTC due to denseness. How to tessellate the cluster to 
achieve a larger MTC will be discussed in the following subsection. In addition, there exists another 
time-division multicast method to multicast in Fig. |4l Namely, the typical transmitter multicasts a 
smaller region in each time slot. We can show that this time-division multicast method attains a less 
Ae than the multihop multicast method. Since it just uses the same transmitter to multicast, the average 
distance from the transmitter to the intended receivers is longer than that of the multihop multicast 




(22) 



Proof: See Appendix III-D[ 



□ 
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method. So time-division multicast has a larger path loss so that the receivers have a lower SIR and 
thus less cluster transmissions are allowed. 



B. Capacity Gain Achieved by Multihop Multicast 

According to Theorem [2l the MTC with multihop multicast obtained from (|22)) and h = 9(log(/c)/r) 



can be concluded as follows 

Since just comparing (|23l) with its single-hop counterpart (l20l) is hard to see if multihop multicast 
achieves a larger MTC or not, defining a capacity gain as follows can help us understand when 
multihop multicast is better than single-hop multicast. 

'a in m\) 



gdv) = 10 log 



10 



(dB). (24) 



C, in mi 

Since the capacity gain is dependent on f , we can formulate an optimization problem with constraints 
on V as follows. 

min —gc{v), subject to 1 — i; < and t> — r < 0. (25) 

V 

Hence, if there exists a minimizer v such that the minimum of —gc is negative, then multihop multicast 
achieves a higher MTC than single-hop multicast. The problem in (l25l) is a convex optimization 
problem as shown in the following theorem. 

Theorem 3. Suppose a packet is delivered by the multihop multicast method and each cluster is 
tessellated into v smaller multicast regions of equal area. (l25l) is a convex optimization problem and 
thus there exists a unique feasible solution of v to it. 



Proof: We have to verify that the optimization problem in (1251) is convex and it has only one 
optimal solution of v. In other words, we have to show that —gc is strictly convex. Now considering 
the noncooperative receiver case, the capacity gain in (|24)) can be explicitly expressed as follows: 



10 



( V^) logio (7) + I logio (1 + 7) + (7 + 2) logio(^) 
z logio (^) + z logio(l + z) + {2 + z) logio(2;) + 2 logio(r) + ^ logio 



ed+r) 
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where z = ^. Taking the derivative of the above equation with respect to z twice and letting < 0, 
it yields the inequality {l + T/v)^y{2r — v)/{2t + v) > 1, which is always true for all r G N+. Thus, 
—gc is strictly convex. □ 
Theorem [3] indicates that the proposed multihop multicast method can achieve a larger MTC than 
single-hop multicast since gc is strictly concave. The reason that multicasting in the smaller regions 
hop by hop can increase the MTC is because the average number of intended receivers in a tessellated 
region and path loss are largely reduced, and thus merely few retransmissions are able to make the 
multicast outage probability lower its designated upper bound. So the increase in is over the loss 
in spectral efficiency. For time-division multicast, no capacity gain is achieved by it since its capacity 
gain function is also negative convex. 

V. Conclusions 

Multicast transmission in a wireless ad hoc network is modeled by a Poisson cluster process, where 
transmitters in the network follow a stationary PPP and each of them is associated with an area-fixed 
cluster in which the intended receivers are a another stationary PPP. The MTC is defined under the 
constraints on the multicast outage probability and the decoding delay. Three network conditions, 
dense, large, and large dense, are specified in order to attain scaling characterizations on MTC. 
The scaling behaviors of the single-hop and multihop MTCs under the three network conditions are 
presented by a general expression. They are affected significantly by the decoding delay but not by 
the fading channel models. In addition, for multihop multicast our main result shows that the MTC 
is superior to its single-hop counter part if all clusters are tessellated appropriately. 

Appendix I 
Useful Propositions 

A. Moment Generating Functional of Stationary Independent PPPs 

Proposition 1. Let i3(0, r) be a circular disc centered at the origin with radius r and = {{Xij, Hij) : 
Xij G B{0, r) nM."^ , r > 1 , j E N} be a stationary marked PPP of intensity Aj for all i E [1, 2, ■ ■ ■ , L] 
and {s/Hij} are i.i.d. Nakagami-m random variables with unit mean and variance. Suppose $j has 
a Poisson shot generating function li : x R_|_ — )■ M_|_ which is defined as li = Xlx e* 
where a > 2. If {^i} are independent, then the sum of the Poisson shot generating functions, i.e. 
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/ = Yli=i ^i' following moment generating functional: 



where 



£,(0i) = E[e-*i^] =exp ^-7rAi(0i, r) 5^ , V01 GM++ (26) 
MM2) = E [e<^^^] = exp l^vr A2(02,r) A,, j , V02 G (0,mr°) . (27) 



Proof: Since the marks {-f^i^} for all points {Xij} in their corresponding PPP are i.i.d. Nakagami- 
m random variables with unit mean and variance, their probability density function is given in Q 
and rewritten in below for convenience: 

1 m 



Thus the Laplace transform of fnih) is 

POO 

^0 



mm roo 



T{m) Jo (1 + w/m)"^' 

Moreover, the Laplace transform for a Poisson shot process $j with i.i.d. marks {Hij} is given by 

m 

CjXh) = exp (-A,, / (1 - E^, [e-<^^^^(l^l)]) dX j , (30) 

I JB(0,r) J 



and we also know 



POO 

Eh [e-*i^^(l^l)] = / e-<^i^(l^l)V//(/i) d/i = Ch (0i^(|X|)) . (31) 
Jo 
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Substituting ([33 into dSO]) and it follows that 



exp -Ai 



B(0,r) 



1- 1 + ^£(|X|) 



exp ( — 27rAj 



1 - 



dX 



exp ^— vrAj 
exp(-7rAi Ai(0i,r)) 



m— 1 



- E 



m 



X dx 



mr°'/4>l p-l+^ 



[i + ty 



■dt 



(32) 



Since all the PPPs are independent, it yields the following desired result: 

1=1 \ i=l J 

Now consider the case of A^/(02)- Similar to Aii^{(p2) can be written as 



A^,,(02)=exp(A. / (l-E^[e^^^^(l^l)])dFl. 

I JB(0,r) ) 



(33) 



We also know Mniw) = J^e^^fHih)dh = j^z^;^ if w E (0,m) and E^,[exp(02i^^(|l^|))] 
■MH{<p2^{\y\)) because 02 G (0,mr"). Then following the same steps in (|32|) . we can show that 



/ / 02 \ " \ ^ / m 



7,(02) = exp ^ttA 

= exp (vrAj A2(02,r)) . 
Therefore, A^/(02) = exp (n A2(02, r) X;f=i K^- 



E 



(34) 

□ 



B. The Duality between PCP and PPP 

Proposition 2. A stationary Poisson cluster process (PCP) can be constructed by a given stationary 
Poisson point process (PPP). Similarly, for a given Poisson cluster process (PCP) a Poisson point 
process can be constructed from it as well. 

Proof: The first statement is based on the definition of a PCP lfT3l . Consider a stationary PPP 
of intensity Ad is given. Then we replace each point Xj G with a random finite set of points 

which is called the cluster associated with point Xj. Then the superposition of all clusters yields 
the stationary PCP Z'^ = IJxje<i><^ ^^j- Suppose now a stationary PCP Z'^ is given. In the following 
we are going to show that a PPP can be constructed by randomly selecting a point in each cluster. 
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Let ^ be a bounded Borel set, be the PPP formed by all parent points in Z^, and <i>'^ consists 
of the points randomly selected from all clusters of Z'^ . The capacity functional of is defined as 
follows [HI: 

Td(^) U(<i>^) > oj , (35) 

where A{^'^) denotes the numbers of points of in a bounded set A C M^. Since the capacity 
functional of a point process completely characterizes its distribution, is a stationary PPP if we 
can show 

%{A) = P[^($^) > 0] = 1 - exp(-Ad/i(^)), 

where fi{A) denotes the Lebesgue measure of set A. 

Consider a large bounded Borel set C such that A C C and C encloses all clusters of Z'^. So C($'^) 
is a Poisson random variable of parameter X^fi^C). Since is formed by the points that randomly 
selected by all points in the probability that all points in ^'^ are not in A is 



P 



A{¥ 







P 







P 







-p 



Ai¥) = 0|C(<l^) 

^(<l^) =0|C(l>^) 7^0 P C{¥)y^O . 

Since is generated from $^ and C covers all clusters of Z^, P C{^^) = = P[C($'^) = 0] and 
thus it follows that 

pU(l>^) = oj = 1 ■p[c($'') = 0] + p U(<l'^) = o|c($'') ^ oj (1 -p[c(<i>^) = 0]) 

= exp(-AdM(C)) + P \Ai¥) = 0|C($^) 7^ Ol (1 - exp(-AdAi(C))). 



Now letting yu(C) — )■ oo (i.e. considering an infinitely large network), we can have 

\,fi{A) Y'"'^'^ 
AdMC)y' 

Therefore, Ta{A) = 1 - exp(-Ad/i(.4)) and $^ is a stationary PPP 



P 



Ai^"^) = 0\Ci¥) ^ 



lim I 1 



exp(-AdAi(^)). 



□ 



Appendix II 
Proofs of Lemmas and Theorems 

A. Proof of Lemma [7] 

The Laplace functional of the stationary PPP $ for a nonnegative function (7 : — )■ 
and shown as follows [fT9l: 



is defined 



CM = E 



■ ^2 g{X) 4.{dX) 



exp 



(l_e-^W)A,MdX) . 
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Since the Laplace functional completely characterizes the distribution of the point process, we can find 
the intensity of by looking for (51). Recall that = {Yj e $[, : maxteii,- ,r] Hj{t)i{\Yj\) > (3Iq} 
and {Hj(t)} are i.i.d. Vt G [1,2, ■■ ■ ,r]. Let be an indicator function which is equal to 1 if 

X E A and 0, otherwise. The Laplace functional of for g{Y) = g(Y)l^c(Y) is given by 

^ i f / (e-3(^)p[F G + 1 - P[r G Ar fiidY)] 

= exp y (1 - e~»(^)) P[F G $'=]Ar ^i{dY)^ . 

Also, for all r G [1, s] and considering At = 6(e), we have 

P[rG$T = F[H^{m\Y\)>f3Io{t),t=l,...,T] 

l-{F[Hi{\Y\)<pio]r 

= l-(E,[Fj,(/3r"Jo|Jo)]r- (36) 

where (a) follows from the fact that the temporal correlation of interference can be neglected for 
small At [|20l . Thus, we have 



CM = exp y' (1 - e-^W) Ac(r, r) /x(dr)^ 



where Ac(r, r) = Ar(l - (E/[FH(/5r"Jo|/o)])^). From the above result we know that $^ C $^ is a 
nonhomogeneous PPP because its intensity Ac(r, r) is the intensity Ar of $q scaled by (l36l) . 
First consider is Rayleigh fading. We can have the following: 

E,[F^,(/3r"Jo)|/o)] = / FH{f3r"co)fj,{u)du 

Jm.+ 

Jo 

where (6) follows from the result of Proposition [T] in Appendix IT A[ So using (|26l) and (l28l) . it follows 
that 

A,(r, r) = Ar {1 - [1 - exp (-7rAi(/3r", oo)At)]^} . 
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Now if yH is a Nakagami-m random variable, then we can have the following: 

E,[F^(/3r%)|/o)] = / Fu{fir'^uj)fj^iuj)d.uj 

Jr+ 

r(m, mPr°'u) 



T(m) 



-//(,(u;)dw, 



where T(y,x) = t^^^e^* dt is the incomplete Gamma function and r(m,0) = T(m) = (m — 1)!. 
Also, if fw{w) is a probability density function of random variable W then we can have the following 
result for m > 1: 



r(m, aw) fw{u!) dw 



/ / t'^-^e-' fwiw)dtdw 

Jo J aw 



(-l)™-ia™ 



im—l 



da""-! V « 

According to Proposition [U in Appendix II-A[ we can have Cig{(f)) with Ai(0, oo). Thus, 



Ej [FH(/3r%|/o))] 



1 + 



r(m) d0™-i 



^/o(0) 



=ml3r°' 



Substituting (|38l) into (1361) . then the result in (fT2)) can be arrived. 



(37) 



(38) 



B. Proof of Theorem [7] 

Here we only provide the proof for Nakagami-m faidn with m > 1 since the proof for Rayleigh 
fading is similar. According to the outage probability ^ upper bounded by e, we know 

^r[\o{R. r)] >\ + ^^^^^ = (l - I) + 6(6^), 

for sufficiently small e. Using the upper bound in (fT5l) and the lower bound on E/j[Ac(-R, r)] obtained 
in above, it follows that 

ER[(l-^('"-i)(m/3i?"))]"<^. (39) 



Therefore, we further have 



Er[^ 



(m-1). 



ml3R')] > 1 



(40) 



Note that Er [^("-^)(m/3 is upper bounded by one and thus it approaches to unity when e/kis 
sufficiently small such that \ftfk < e. That means \E''^™~^)(m/3 -R") is very close to 1 almost surely and 
thus At = 0(e). If k is sufficiently large, then we have exp(— tt AtAi) = 1— tt AtAi+0(e^). Substituting 
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this expression into ^(™-i)(m/3i?") can be reduced to 1 - TrAt Ai njLi^(l - 2/ja) + 0(e 



Define Ai(0) = [Ai(0, oo) - Ai(0, a)] lYj=ii^ - 2/jtt)- Choosing < a < 1, we have 



(r + 1] 



According to (1391) and the above result, these two upper bounds should coincide when the maximum 
contention intensity At is reached. Namely, (fT6l) is obtained from the following: 

^e(r + l) 



A, 



7rs2/3l v^Ai(/3) 



+ e(e^) ^ A, 



r + 1 



us- 



73^ Vk' 



(41) 



where r] = l/Ai(/3). This completes the proof. 



C. Proof of Lemma \3\ 

According to Campbell's theorem, we know E[/o] = and then the lower bound on h can be 



reduced as follows. 



(a) 

6 > E 



log 1 + 



H-n 



(b) 

> E 



log 1 + 



H„ 



(«-2) 



27rs2A, 



> log 1 + 



TTS^At 



At / \a — 2^ 

where (a) follows from Jensen's inequality since log(l + a/x) is convex for a; > and constant 
a > 0, (6) is due to s"^^ > 1, and Ff^^^^ is the CCDF of random variable -f^max- In addition, we also 
know h in (ITTl ) is upper bounded as follows. 

h < log (E[/o] + E[i7^ Js-") + E [log (l//o)] . 

The upper bound is obtained by first conditioning on Jq and then using Jensen's inequality. We know 
the term log(l//o) in the upper bound is convex and thus E[log(l//o)] > log(l/E[/o]). Since the 
network is interference-limited, E[/o] must be bounded above zero. Thus there exists a 71 > such that 
log(7i) < E[log(/o)] < log(E[/o]), which means E[log(/o)] > log(72E[/o]) for any 72 G (0, 7i/E[Jo]]. 
So the upper bound can be simplified as follows: 



log (E[/o] + E[/J^ax]s-") +E 



log I — 

-'0 



< l0£ 



E[/o] 

1 



+ E[i?- 



log 



72E[/o] 
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because s"" < E[i?-"] < ^2(^ and E[H^^^] > 1. Note that the CDF of H^^^, Fh^,,, is equal to 
{FffY since {Hj{t),t = 1, . . . , r} are i.i.d. random variables. The proof is complete. 

D. Proof of Theorem |2] 

We start with the derivation of the upper bound on E[Ac] as follows, (i) The upper bound on 

E[Ac]: Suppose the multicast cluster IZq is tessellated into v smaller multicast regions of equal area 
-Ks^ /v. Let T denote the set of all the pathes from the source multicast region (i.e., the tessellated 
region in which the central transmitter is located) to the last multicast region. Let T* G T be the 
chosen path for multicasting a packet and Si denote the iih. tessellated region on T^. Since there are 
V multicast regions on T*, the duration of a time slot for multicasting a packet in each region is at 
most T /v transmission attempts. In order to implement multihop multicast, the multicast outage event 
should be redefined here as follows. We say there exists a multihop multicast outage in a multicast 
time slot if any of the receivers in the current multicast region does not receive the packet or no 
receivers in the next multicast region receive the packet after t/v transmission attempts. In addition, 
is chosen according to the following criteria: 



where Si is the multicast outage event happening in Si. So is the path with the minimum probability 
of end-to-end multicast outage. 

If is the receiver-connected process in Si after r attempts, then Si can be expressed as C 
($,.n5i)U(yi+i G Si^i)} where Kj+i is any intended receiver in and it could be the transmitter in 
the {i + l)th time slot. Due to the broadcast nature of wireless communication, the intended receivers 
in the regions other than Si could also receive the packet which is merely multicasted to Si in the 
zth time slot. Thus, if the average intensity of $s- at the end of the tth time slot is Ex[As-(X, tr/t>)] 
where X G 7^o, then we know Ex[A,,(X, tr/t;)] < Ex[A,2(X, tr/t;)] < ■■■ < Ex[A,„(X, tr/t;)] 
since is a filtration process. 

2 

Since the Lebesgue measure of 5j is — , the probability of Ei on is given by 



V 



= arermin 

Ter 



i-ll{i-m\T]) 



(42) 



i=l 
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where = 1 — (1 — e)^/^ is the upper bound constraint on the multicast outage probability in each 
tessellated region. Considering sufficiently large k, the above equation can be simplified as 

F[£,\T,] = 1 - exp ^ (Ar - A,^(X, zr/v)) fi{dX)^ + Q{e^). 

The multicast outage event in TZq is the union of all the multihop multicast outage events, i.e., 
£ = |jr=i ~ (nr=i ^i^)'^- since all multihop multicast outage events are independent, the probability 
of £ can be explicitly expressed as 

F[£\T,] = 1-Y[F[£^\T,] 



i=l 



= l-exp|-7rs'Ar + 2^y XsXX,tT/v) i2{dX) } + e{e') < e. 
According to ([8]) and the above equation, it follows that 

V ^ — ' /c TTS /V V ^ — ' 

i=l ' i=l 

From the results in Lemma [21 we can know Ex[As.(X, ir/v)] with Rayleigh and Nakagami-m fading 
at the ith time slot. Let A{r) = 1 — exp(— 7rAi(/3r°', oo)At) for Rayleigh fading or A{r) = 1 — 
^(m-i)^^^^a-j for Nakagami-m fading with m > 1, and thus the upper bound on E/j[Ac(-R, r)] can 
be shown as 



(ii) The scaling law of A^: Using the upper bound on E,ji[Xc{R,T)] obtained From ([8]) and the 

above result, it yields the following inequality: 

v^^ r./T.M- Er\A\v(1 -ErIAV) 
k ^ l-Eji[A]v 

Also, we know = 1 - </T^ = ^ + e{e^) and En[A]^ < ^ = f + ©(e^)- Thus, we know 

E^[^('"-i)(m/3i?")] ^ 1 since Efi[^(™-i)(m/3i?")] < 1 and A; > So the random variable 

is very close to one almost surely and thus At is 0(e). Then following the same steps 
in the proof of Theorem [H we can show that 

Thus, (|22l) is arrived, which completes the proof. 



28 



References 

[1] p. Gupta and P. R. Kumar, "The capacity of wireless networks," IEEE Trans. Inf. Theory, vol. 46, no. 2, pp. 388^04, Mar. 2000. 
[2] S. Shakkottai, X. Liu, and R. Srikant, "The multicast capacity of large multihop wireless networks," in Proc. ACM International 

Symposium on Mobile Ad Hoc Networking and Computing (Mobihoc'07), Sep. 2007. 
[3] U. Niesen, P. Gupta, and D. Shah, "The unicast and multicast capacity regions of large wireless networks," submitted to IEEE 

Trans, on Inf. Theory, May 2009. 
[4] X.-Y. Li, "Multicast capacity of wireless ad hoc networks," IEEE/ACM Trans. Netw., vol. 17, no. 3, pp. 950-961, Jun. 2009. 
[5] A. Keshavarz-Haddad, V. Ribeiro, and R. Riedi, "The multicast capacity of large multihop wireless networks," in Proc. Annual 

International Conf. on Mobile Computing and Networking, Sep. 2006. 
[6] B. Tavil, "Broadcast capacity of wireless networks," IEEE Commun. Lett., vol. 10, no. 2, pp. 68-69, Feb. 2006. 
[7] P. Chaporkar and S. Sarkar, "Wireless multicast: theory and approaches," IEEE Trans. Inf. Theory, vol. 51, no. 6, pp. 1954-1972, 

Jun. 2005. 

[8] P. Jacquet and G. Rodolakis, "Multicast scaling properties in massively dense ad hoc networks," in Proc. International Conf on 

Parallel and Distributed Systems, Jun. 2005. 
[9] A. Keshavarz-Haddad, V. Ribeiro, and R. Riedi, "Broadcast capacity in wireless multihop networks," in Proc. ACM International 
Symposium on Mobile Ad Hoc Networking and Computing (Mohihoc'06), Sep. 2006. 
[10] R. Zheng, "Information dissemination in power-constrained wireless network," in Proc. IEEE INFOCOM, Apr. 2006. 
[11] A. Keshavarz-Haddad and R. Riedi, "On the broadcast capacity of multihop wireless networks: Interplay of power, density and 
interference," in Proc. Annual IEEE Communications Society Conf. on Sensor, Mesh and Ad Hoc Communications and Networks 
(SECON'07), Jun. 2007. 

[12] S. Weber, X. Yang, J. G. Andrews, and G. de Veciana, "Transmission capacity of wireless ad hoc networks with outage constraints," 

IEEE Trans. Inf. Theory, vol. 51, no. 12, pp. 4091^102, Dec. 2005. 
[13] D. Stoyan, W. Kendall, and J. Mecke, Stochastic Geometry and Its Applications, 2nd ed. John Wiley and Sons, Inc., 1996. 
[14] R Baccelli, B. Blaszczyszyn, and P. Miihlethaler, "An Aloha protocol for multihop mobile wireless networks," IEEE Trans. Inf. 

Theory, vol. 52, no. 2, pp. 421-436, Feb. 2006. 
[15] M. Haenggi, J. G. Andrews, F. Baccelli, O. Dousse, and M. Franceschetti, "Stochastic geometry and random graphs for the 

analysis and design of wireless networks," IEEE J. Sel. Areas Commun., vol. 27, pp. 1029-1046, Sep. 2009. 
[16] O. Arpacioglu and Z. J. Haas, "On the scalability and capacity of wireless networks with omnidirectional antennas," in Proc. 

Information Processing in Sensor Networks (IPSN'04), Apr. 2004. 
[17] O. Dousse and P. Thiran, "Connectivity vs capacity in dense ad hoc networks," in Proc. IEEE INFOCOM, Apr. 2004. 
[18] J. F. C. Kingman, Poisson Processes. Oxford University Press, 1993. 

[19] F. Baccelli and B. Blaszczyszyn, Stochastic Geometry and Wireless Networks. Now Publishers: Foundations and Trends in 
Networking, 2009. 

[20] R. K. Ganti and M. Haenggi, "Spatial and temporal correlation of the interference in Aloha ad hoc networks," IEEE Commun. 
Lett, vol. 13, no. 9, pp. 631-633, Sep. 2009. 



29 



Intended * • \ 

Mul^cast -n' 

Area a / 

(cliister) • / y |, 

V —X L- 



/ RX A ^ ' ' « 



\ 

• > 



Fig. 1. Multicast transmission model in a planar wireless ad hoc network: The transmitters in the network form a stationary PPP 
of intensity At. Each transmitter (triangle) has an intended multicast area of radius s, where all the intended receivers (small circles) 
also form a stationary PPP of intensity Ar. A transmitter and its corresponding intended receivers are indicated by the same color in 
a cluster. So each cluster could contain other transmitters and unintended receivers in addition to its own transmitter and intended 
receivers. 
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Fig. 2. The simulated MTC of a large network with Rayleigh fading for e = 0.1, /? = 2, a = 3 and Ar = 0.1. 
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Fig. 3. The simulated MTC of a large dense network with Rayleigh fading for e = 0.1, /3 = 2 and a = 3. 




Fig. 4. Each multicast cluster in a planar network is tesselated into 6 smaller multicast regions of equal area. The arrows in each 
cluster show an example path of delivering a packet over the tessellated regions by the multihop multicast method. 
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Fig. 5. The simulated MTC achieved by multihop multicast in a large network with Rayleigh fading for e = 0.1, /3 = 2, q = 3, 
T = 20 and Ar = 0.2. 
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Fig. 6. The simulated MTC achieved by multihop multicast in a large dense network with Rayleigh fading for e = 0.1, /3 = 2, a = 3 
and T = 20. 



